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Abstract. We show that the Circular Orthogonal Ensemble of random matrices arises 
naturally from a family of random polynomials. This sheds light on the appearance of 
. . . random matrix statistics in the zeros of the Riemann zeta-function. 
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^ ■ 1. Introduction 

! The statistics of eigenvalues of unitary matrices, chosen uniformly with respect to Haar 

measure on U(N), are observed to closely match the statistics of zeros of the Riemann zeta- 
function and other L-functions |2~T1 122 EU ■ in addition, eigenvalues of other compact 
classical matrix groups give a good model of the zeros of various families of L-functions |161 
! H2J I2H EH IE] ■ Furthermore, the characteristic polynomials of the matrices provide a good 

model of the L-functions themselves El EH HI EH] • 

The L-functions studied in number theory are Dirichlet series having a functional equation 
and an Euler product. In this paper we are concerned with a wider class of Dirichlet series 
which have a functional equation but do not have an Euler product (see, for example, [9, 
EH])- It has been suggested that such functions can be modelled by random self-reciprocal 
polynomials |5j. In Appendix A we discuss the example of Epstein zeta functions. It is 
Q^. also possible to create Dirichlet series with a functional equation but no Euler product from 

Maass forms on a nonarithmetic group [7j. Since these functions do not have an Euler 
product, they are not expected to satisfy the Riemann hypothesis. However, it is possible 
that occasionally such functions will satisfy the Riemann hypotheses or, as is more likely, 
will have a large number of consecutive zeros, say the first 100 trillion of them, on the critical 
line. 

This paper is motivated by the following questions: do the Riemann zeta-function and 
the other L-functions of number theory behave differently than random Dirichlet series with 
functional equation which just happen to satisfy the Riemann hypothesis? That is, does the 
Euler product have any effect on the zeros beyond forcing them to be on the critical line? 
Our results suggest that the answer is 'y es \ and the Euler product also has an effect on the 
local statistics of the zeros. 

1.1. Random polynomials. In random matrix theory, the characteristic polynomial of a 
random matrix can be viewed as a random polynomial where the randomness is explicitly 
encoded in the zeros. For example, consider the Weyl integration formula for the classical 
compact groups [2E1E3 or the j3 = 1, 2, 4 ensembles of random matrix theory [20J. On the 



Date: February 7, 2008. 

Research supported by the American Institute of Mathematics and the Focused Research Group grant 
(0244660) from the NSF. This work was started during the program on Random Matrix Applications in 
Number Theory at the Isaac Newton Institute for Mathematical Sciences. The second and third authors 
were also supported by a Royal Society Dorothy Hodgkin Fellowship, and the third author was partially 
supported by EPSRC. 

1 



2 DAVID W FARMER, FRANCESCO MEZZADRI, AND NINA C SNAITH 

other hand, usually in the study of random polynomials (for a review of the subject see, e.g., 
Farahmand the randomness is explicitly encoded in the coefficients of the polynomials. 
We will show that there is a simple but surprising connection between these two perspectives, 
and we believe this connection is relevant to the appearance of random matrix statistics in 
the zeros of the Riemann zeta-function and other L-functions. 
The polynomials we consider are of the form 

(1.1) f(z) = z N + a lZ N - l + --- + a N 
with |aj\r| = 1 which have the symmetry 

(1.2) f(z) = a N z N J 

where f(z) := f(z). Such polynomials are called self -reciprocal. Equation ()1.2|) ensures that 
the zeros of / occur either on the unit circle or in pairs located symmetrically with respect 
to the unit circle. The symmetry ()1.2|) is analogous to the functional equation of L-functions 
which arise in number theory, see ()1.12j) below. In terms of the coefficients cij the functional 
equation (jl.2|) translates into the relation 

(1.3) CL]\f—j = CL]\[Clj. 

The change of variables z = e lx in the function a N 2 z~ N ^ 2 f(z) converts a self-reciprocal 
algebraic polynomial into a real trigonometric polynomial 

(1.4) c n cos({N/2-n)x) + d n sm((N/2-n)x) , 

0<n<N/2 

where c n and d n are real. And if / has real coefficients then the associated trigonometric 
polynomial is an even function of x, having only cosine terms in its expansion. The roots of 
these trigonometric polynomials are either on the real line or in complex conjugate pairs. This 
is somewhat closer to the symmetry of an L-function, but in this paper we phrase everything 
in terms of algebraic polynomials (J1.1J) to emphasize the comparison with characteristic 
polynomials of random matrices. 

For even trigonometric polynomials, whose coefficients are independent standard normal 
random variables, Dunnage |2j discovered that the expected number of real zeros is given by 




(1.5) ^ + O(A^/ 13 (logA0 3 / 13 ). 

v3 

Bogomolny et al. [2] extended Dunnage's result and studied the average fraction of roots of 
the polynomial (jl.lj) on the unit circle when = 1 and the coefficients ctj (for 1 < j < N/2) 
are complex normal random variables with mean zero and whose standard deviation a varies 
with N as 

(1.6) o = e/y/N. 

They discovered that in the limit e — > 0, i.e. when all the coefficients have a narrow 
distribution centered around zero, there is a small neighbourhood of e where the average 
fraction of roots lying on the unit circle is one. Instead, when the coefficients have a broad 
distribution, i.e. in the limit e — > oo, the average fraction of zeros on the unit circle tends 
to l/y/3. They also computed the two-point correlation function of the zeros on the unit 
circle, which for short distances grows linearly. 
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1.2. L-functions. An L-function is a Dirichlet series, 

oo 
n=l 

with a n = O t (n e ) for every e > 0, which has an analytic continuation to the complex plane 
(except for possible poles on the line Res = 1) along with two additional properties. First, 
it has a functional equation 

(1.8) := lL (s)L(s) =e5(l-s), 
with |e| = 1 and of the form 

w 

(1.9) lL (s)=P(s)Q s l[r(w,s + fi ] ), 

3=1 

where Q > 0, Wj > 0, Re/ij > 0, and P is a polynomial whose only zeros in a > are at the 
poles of L(s). Second, it has an Euler product representation of the form 

(1.10) L(s) = l[L p (l/p s ), 

v 

where the product is over the primes p, and 

oo 

k=0 P 

where b n = 0{n e ) with 9 < |. It is conjectured that such L-functions satisfy the Riemann 
hypothesis, which is the assertion that the nontrivial zeros lie on the critical line Re(s) = |. 

The functional equation is a symmetry with respect to the line Re(s) = |, so the nontrivial 
zeros of L are on the line Re(s) = ~, or they are located symmetrically on either side of it. 
To show the analogy with self-reciprocal polynomials ()1.2|) it is more convenient to write the 
functional equation in asymmetric form: 

(1.12) L(s)=eX L (s)L(l-s), 

where X L (s) = 7^(1 — s)/7i(s). Note that \X L \ = 1 on the 1-line, that is, on the line of 
symmetry of the L-function. Thus, there is a perfect analogy between the self- reciprocal 
property of / and the functional equation of L. 

1.3. Results and discussion. For Dirichlet series with a functional equation but no Euler 
product, it is not expected that all zeros lie on the critical line. However, it is possible that 
even if such a Dirichlet series does not satisfy the Riemann hypothesis, it could have, say, its 
first 100 trillion zeros on the critical line. Our motivating question is this: is the Riemann 
zeta-function, with its Euler product, distinguishable from a random Dirichlet series with 
functional equation that just happens to have its first 100 trillion zeros on the critical line? 
By considering the analogous case of random self-reciprocal polynomials, we suggest that 
the answer is 'yes'. 

Let us consider the space IZn of all the self-reciprocal polynomials of degree N . Given a 
nonvanishing probability distribution onTZ^ (meaning that nonempty open sets have positive 
measure) the subset Cn of those polynomials which happen to have all of their zeros on the 
unit circle has positive measure. The proof is quite simple. The zeros of a self- reciprocal 
polynomial lie either on the unit circle or in pairs symmetric with respect to the unit circle, 
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and they are continuous functions of the coefficients. Two zeros on the circle must "collide" 
in order to move off the circle, so there is a small open neighborhood of coefficients in which 
the zeros remain on the circle. This neighbourhood has positive measure. Therefore, the 
restriction to of the measure on 1Z N is unique and well defined. Such a restriction is then 
made into a probability measure by renormalizing the volume of Cat to one. 

Our interest is mainly concentrated on the distribution of the roots of those polynomials 
whose zeros lie all on the unit circle. Therefore, our approach will be to define a joint 
probability density function on the coefficients of the polynomials in TZ^ and to study its 
restriction to Cjv- The coefficients a n of polynomials with all their zeros on the circle are 
bounded in magnitude by (^), therefore not only has Cn positive measure, but it is also 
compact. Thus, the most natural choice is to put a distribution on the coefficients a n 
which is uniform on a bounded disk containing Cjy and zero outside. The following theorem 
gives the joint probability density function for the roots of such polynomials. Petersen and 
Sinclair [23J have found some interesting geometric properties of the coefficients of these 
polynomials. They have also proved independently a result very similar to Theorem II .11 In 
their work (Lemma 4.1 of [23J) they fix the coefficient a at, rather than allowing it to vary, 
but they arrive at the same Vandermonde form for the Jacobian. 

In what follows we denote by A(xi, . . . ,xn) the Vandermonde determinant, i.e. 

(1.13) A(x 1 ,...,x N ) = Y[i x k -Xj), 

j<k 

and we denote by e n the nth elementary symmetric function 

(1.14) e n (x 1 ,...,x N )= ^2 x h'-- x in- 

l<ii<— <i n <N 

Theorem 1.1. Suppose N is odd. Consider random monic polynomials z N + Y2n=i a n zN ~ n 
satisfying the self-reciprocal property a^_ n = aN~a n , with di, . . . , 0(at-i)/2 chosen indepen- 
dently and uniformly in \a n \ < ( ) ; and with <p chosen uniformly in [0, 27r), where On = ^ , 
and restrict to those polynomials having all zeros on the unit circle. The joint probability 
density function of the set of zeros e tSl , . . . , e l5N is given, up to a normalization constant, by 

(1.15) \A{e i5 \...,e i5N )\ = JJ \e i5k -e^|. 

j<k 

For N even, consider self-reciprocal random monic polynomials with a\, . . . , ajv/2 chosen 
independently and uniformly in \a n \ < CT), arid restrict to those polynomials having all zeros 
on the unit circle. The joint probability density function for the set of zeros e lSl , . . . , e l&N is 
given, up to a normalization constant, by 

(1.16) \e N/2 (e l5 \ . . . , e"*)A(e*V . . , e*»)| = \e N/2 (e i5 \ . . . , e*»)| J] |e** - e*>\, 

j<k 

where e N/2 {e i5 \...,e i5N ) is the N/2 elementary symmetric function in the variables e l , . . . , e 
and is equal to (—1) N ^ 2 times the coefficient a^/ 2 . 

In particular, the joint probability density function for odd N is the same as that for 
eigenvalues of a randomly chosen matrix in the Circular Orthogonal Ensemble COE(N). 

Some aspects of this result can be derived from |2Zj but they approach the subject from 
a different perspective. 
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The theorem suggests that if a Dirichlet series with functional equation is chosen at ran- 
dom, and all of the zeros in a particular interval of the critical strip happen to lie on the 
critical line, then those zeros should have similar statistics to those of eigenvalues of matri- 
ces from the COE. The COE is the symmetric space U(N)/0(N) with the measure induced 
from Haar measure on U(N). Thus the joint probability density function for the eigenvalues 
is 

(1 17) \A(e iSl e i5N )\ = - TT \e iSk - e iSj \ 

K ' (4V^)*T(l + JV/2) 11 n (4^)^(1 + ^/2) 11 16 

j<k 

All numerical calculations of zeros of L-functions having an Euler product show the statistics 
of the CUE (the group U(N) with Haar measure), so this suggests that the L-functions from 
number theory are not typical Dirichlet series with a functional equation. 

Thus, the Euler product which is considered a necessary condition for a Dirichlet series with 
functional equation to satisfy the Riemann hypothesis does more then just force the zeros 
onto the critical line. The Euler product fundamentally changes the nature of the spacing 
of the zeros, in particular changing the linear repulsion of zeros of random polynomials and 
the COE into the quadratic repulsion of the CUE, whose joint probability density function 
is given by 

i i 

(I 18) \A(e iSl e i&N )\ 2 = TT \e iSk - e iSj \ 2 

In random matrix theory it has long been conjectured that in the limit iV — > oo the 
local correlations of the eigenvalues of random matrices depend exclusively on the invariance 
properties of the probability distribution that defines the ensemble and not on the explicit 
form of the measure itself. This random matrix hypothesis is one of the most important 
features of the subject. Mathematically it translates into the statement that, provided that 
the local eigenvalue density has the same asymptotic behaviour, the local correlations are 
mainly determined by the absolute value of (powers) of the Vandermonde, whose origin is 
essentially geometrical. Therefore, we conjecture that in the limit iV — > oo the local statistics 
of the roots of the polynomials in Cn will be independent of our choice of the joint probability 
density function for the coefficients of the polynomials in TZn- Indeed, because of this reason 
the extra factor |ejv/2 1 appearing in equation (|1.16|) of theorem 11.11 when N is even should 
not affect the local correlations of the roots in the limit iV — > oo. 

Bogomolny et al. j2] studied self- reciprocal polynomials whose coefficients are independent 
complex normal random variables and computed the two-point correlation function R2(S) of 
the subset of zeros that lie on the unit circle. As 5 — > they observed linear repulsion between 
the arguments of the zeros. Such level repulsion is a direct consequence of the Vandermonde 
that appears in equations ()1.15|) and ()1.16|) of theorem 11.11 (or more appropriately for their 
case in equations (|3.1J1 and (|3.2j) of lemma EH} and supports our conjecture. 

The appearance of CUE statistics for arithmetic L-functions has been compared to the 
appearance of CUE statistics in a chaotic system without time-reversal symmetry [T71 . 
Indeed, the appearance of the CUE statistics for zeros of L-functions has been heuristically 
explained by the analogy between the periodic orbit sum for the density of states of a 
classically chaotic system with no time-reversal symmetry and the density of zeros of the 
Riemann zeta function expressed as a sum over primes. Our observation on the effect of the 
Euler product can be viewed as further evidence for that point of view. 
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In their article Bogomolny et al. [2] commented that the linear behaviour of R2{8) in 
the limit 5 — > was surprising, because it is typical of quantum mechanical systems whose 
dynamics is invariant under time reversal. Instead, they would have expected quadratic 
repulsion between the roots lying on the circle, which is characterizes chaotic systems whose 
dynamics is not time reversal invariant and that self-reciprocal polynomials with arbitrary 
complex coefficients were expected to model. In the case that we study we are facing the 
same paradox. We are modelling L-functions having a functional equation but no Euler 
product and observe that the zeros lying on the critical line should have the same local 
correlations as the eigenvalues of matrices in the COE ensemble. Instead, the zeros on the 
critical line of L-functions that have the further constraint of being expressed in terms of 
an Euler product are correlated like the eigenvalues of matrices in the CUE ensemble. The 
natural expectation would be the opposite, since matrices in the CUE ensemble are not in 
any way restricted, except that the measure on their set should be invariant under group 
multiplication. The COE ensemble is obtained by imposing extra symmetry constraints on 
the matrices in the CUE, whose set, therefore, includes the set of matrices forming the COE. 

We also consider the case of real self-reciprocal polynomials. That is, polynomials of the 
form (jl.lj) satisfying (jl.2|) . where the a n are real. These polynomials have their zeros in 
complex conjugate pairs, and in particular the zeros near z = 1 would be expected to have 
somewhat different behavior than the bulk of the zeros. Anomalous behavior of the low-lying 
zeros occurs for families of arithmetic L-functions, and this behavior has been modelled by 
the low-lying eigenvalues of matrices from the classical Symplectic and Orthogonal groups. 

Theorem 1.2. Let N = 2M. Consider random real raonic polynomials z N + Yl n =x a n zN ~ n 
satisfying the self-reciprocal property a n = aN- n , with a n chosen independently and uniformly 
in the interval \ a n \ < f ) , and restrict to those polynomials having all zeros on the unit circle. 
The joint probability density function of the set of zeros e ltl , e~ ltl , . . . , e ltM , e~ ltM is given, up 
to a normalization constant, by 

(i.i9) yi \ eitm - e ~ itm \ n \ eih ~ e * <j i \ eitk ~ e ~ it3 \- 

In the case where the polynomial has odd degree N = 2M + 1 with roots —1, e ltl , e~ ltl , . . . , 
e ltM , e~ ltM then the joint probability density function of e ltl , e~ ltl , . . . , e ttM , e~ ltM is the same 
as given above. 

The above measure can be written as, up to a normalization constant, 
(1.20) 

We note that this is the square root of Haar measure of USp(N). Thus, random real self- 
reciprocal polynomials, restricted to have all their zeros on the unit circle, do show anomalous 
spacings in their low lying zeros. But it is not the same anomalous spacing that has previously 
been found in families of arithmetic L-functions. This suggests that for the low-lying zeros of 
a family of L-functions with real coefficients, the Euler product has an effect on the vertical 
spacing of the zeros, and those L-functions behave differently than random real Dirichlet 
series with functional equation which just happen to have their first few zeros on the critical 
line. 
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In the next section we prove the Theorems and in the following sections we give the 
Jacobian calculations required in the proofs. We thank Christopher Sinclair for helpful 
information. 



2. Proofs of the theorems 

We are given a measure on a set of polynomials described in terms of the coefficients of 
the polynomial, and we wish to describe the measure in terms of the roots of the polynomial. 
Therefore we must compute the Jacobian of the change of variables from the coefficients to 
the roots. It is well known that the Jacobian is the Vandermonde in the roots in the case 
of polynomials with real coefficients, and the Vandermonde squared in the case of complex 
coefficients. Thus we expect the Jacobian to be close to a Vandermonde in the case of self- 
reciprocal polynomials, but we were unable to find all the results we needed in the literature, 
so we give a complete proof below. 

If X is a random variable then we let (X) denote the expected value of X. In our case 
/ will be a random polynomial and X = M[f] is some function of /, and we will need to 
compute (M[f]). 

In Theorem 11.11 consider N even so a±, . . . , a N / 2 determine /. If p C oeffs( a ij • • • > a N/2) is a 
probability measure on the coefficients of /, which is supported on the set S, then 

(2.1) (M[f]} = J M[/]p coeffs (ai, . . .,a N / 2 )dai ■ ■■da N/2 , 

s 

and by definition f p C oefis(oij • • • , dN/2)da\ ' ' ' da N / 2 = 1. 
s 

We can express (M[f]) in terms of the zeros z\, . . . , Zn of /. Set 

(2.2) Pzeros^l, ■ ■ - ,Z N ) = p cocSs (til (z\ , . . . , 2jv), • • • , a N {z X , • • • , Z N )) 

and let Jc{z\, ■ ■ ■ , Zn) be the Jacobian of the transformation from the coefficients to the 
zeros. Then 

(2.3) (M[f]) = J M[f]p zems (z 1 ,...,z N )\J c (z 1 ,...,z N )\dz 1 ---dz N , 

S' 

where 5" is the image of S under the coordinate change. To prove the Theorems we merely 
specialize this discussion to our particular cases. 

Proof of Theorem } 1. 11 We have that p C oeffs is constant and S is the set of coefficients of self- 
reciprocal polynomials having all their zeros on the unit circle. So p zer0 s is also constant and 
S' = S 1 x • • • x S 1 where S 1 is the unit circle. Thus, we only require the Jacobian of the 
transformation, which is given in the following Lemma. 

Lemma 2.1. If the self-reciprocal polynomial f{z) = z N + ^2n=i a n zN ~ n ^ as a ^ ^ s zeros 
on the unit circle, e lSl , . . . ,e t5N , then when N is odd, the absolute value of the Jacobian of 
the transformation from coefficient variables Reai, Imai, . . . , Rea^N-i)/2, Ima(jv-i)/2, (where 
Ojv = ^ ) to the zero variables 5i, . . . ,5n is given by 



(2.4) 
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When N is even, the absolute value of the Jacobian of the transformation from coefficients 
Reai, Imai, . . . , Rea N / 2 , Ima7v/2 to zeros 5±, . . . , 5n is given by 

(2.5) |j c (e^,. . . , e^)\ = 2~f \e N/2 (e l5 \ e^) A(e l5 \ . . . , e^)\ , 

where e N/2 {e i5 \...,e i5N ) is the N/2 elementary symmetric function. 

A proof of the Lemma can be found in Sectional 
Assembling the pieces we have, for example for odd N, 

■y p2tt p2tt 



(2 - 6) mf] >=j4vwwTW)i -y. ^i^.-.^i*-^. 

which is equivalent to Theorem ll.il □ 

The proof of Theorem 11.21 is identical except that we require the following lemma, which 
is proven in Section HJ 

Lemma 2.2. If the real self -reciprocal polynomial f has even degree N = 2M and has all its 
zeros on the unit circle, then the absolute value of the Jacobian of the transformation from 
coefficients a\, . . . , om to zeros e ltl , e~ ltl , . . . , e l<M , e~ ltM is given by 



[2.7) \J R {e ltl ,e- Ul ,...,e itM ,e- ltM ' 



Y[(e itm - e~ itm ) Y[(e itk - e^)(e itk - e~ u ^ 

m j<k 

In the case the degree N = 2M + 1 is odd, and f has zeros at —1, e ltl , e~ ltl , . . . , e ttM , e~ ltM , 
the Jacobian is again given by the above formula. 

3. Calculation of the Jacobian: complex case 

We prove the following generalization of Lemma 12.11 

Lemma 3.1. Let the roots of a self-reciprocal polynomial f be a\ — e tSl , . . . ,az = e l&L , for 
those roots on the unit circle, and Pi = pie 101 , =- = . . . , Pm = Pm^ 9m , j— = ^j- for the 
roots occurring in pairs off the unit circle. When N = L + 2M is odd, the absolute value of 
the Jacobian of the transformation from coefficients Reai, Imai, . . . , Rea(jv-i)/2j ^- m o J (N~i)/2, <fi 
(where a N = e 1 ^ ) to zeros pi, 9i, p 2 , 2 , . . . , p M , 9m, o~i, ■ ■ ■ ,$l is given by 

(3.1) | Jc(pi,9i,P2,02, ■ ■ ■ , Pm,0m,5i, ...,5 l )\ 



2Af-(iV-l)/2 



i Pm-} 15 5 ^1) • • • 5 

Pm 



When N = L + 2M is even, the absolute value of the Jacobian of the transformation 
from coefficients Reai, Imai, • • • , Rea^/ 2 , ImaAr/2 to zeros p%, 9i, p 2 , 9 2 , . . . , pm, 9m, o~i, ■ ■ ■ ,8l 
is given by 

\Jc(Pi,&i,p2,9 2 , . . . ,pm,9m,Si, . . . ,5 l )\ 



(3.2) = 2 M " 



5 Pm-} -75 •} ^1) • • • •} &l) 

Pm 



x A(/9 X , p 2 , ■ ■ ■ , Pm, j^, <*i, • • • , a L ) 
where the ejv/2 is the N/2 elementary symmetric function. 
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Proof. The polynomial f(z) has order N = L + 2M: 



(3.3) 



f(z) = (z - a x )(z -a 2 )---(z- a L )(z - P ± )(z - J-) • • • (z - (3 M )(z - J-) 

Pi PM 



z N + aiz N 1 + a 2 z N 2 + ■ ■ ■ + CLn-2 z2 + &N-\ Z 



Note that if we define 
(3.4a) aj 
(3.4b) Pi 
1 

% 



(3.4c) 



e l0 \ 

Pi 



j = l,...,L 
j = l,...,M 



_1 

Pi 



J 



M 



(with 5j, pj and 9j real) then is on the unit circle and 



a N = (-ire 



N o 2i6 le 2i0 2 



(3.5) 
Also 

(3.6) aj = (-l) 3 e j (a 1 ,...,a L ,(3 1 ,±,...,(3 M ,j^), 

where ej(xi, . . . ,x n ) is the jth elementary symmetric function. In addition we have 

(3.7) Uj^—j = Ojvflj, 

since by construction / is self-reciprocal. 

For now we take iV odd; the slight variation when N is even is described at the end 
of this section. We want the Jacobian of the transformation from the independent real 
variables Reai, Imai, . . . , Reaiv^i , Imajv-i , <p to the real independent variables p±, Q\, p 2 , 9 2 , 

. . . , Pm, ®Mi • • • 5 



So, the Jacobian is: 

<9Reai 

dpi 
dlmai 

dpi 

./• = dRea jv_ 



dpi 
dlma jv-i 



dRcai 
dOi 

dlmai 
d8i 



dRca jv_i 

Wi 
dlma jv-i 



dpi 
dpi 



d0i 



d0i 



dReai 

almai 
dpM 



dRea jv-i 

dp M 
dlma jv-1 



dPM 
8<t> 

dp M 



9Rcai 

96 M 
plmai 

de M 



dRca jv-i 

„ dd M 
dlma jv-i 

9^ 



(3.8) 



= H) 



JV-1 
2 



dai 
dpi 

dai 
dpi 



dai 
d6i 

dai 
ddi 





96 M 


dai 


dai 


dpM 


dd M 


dai 


dai 


dpM 


90m 



da i 



da i 



da. 



da i 



dReai 

dSi 
dlmai 

d&i 



dRca jv-i 

Wi 

Olmajv-i 

d5[ 
d<t> 
dSi 
dai 
dSi 

dai 
dSi 



da i 



dpi 

da N-i 
dpi 

9<j> 
dpi 



d8i 

da~ N-i 
d9i 

d<t> 
96i 



dpM 

dpM 

d<t> 
dpM 



d6 M 

da. N-i 

de M 

d<f> 

de M 



dSi 

da N-i 
d&i 

dSi 



8Rcai 
d6 L 

dlmai 
d5 L 



dRca i 



dS L 
9Imajv-i 



dS L 
d<f> 
dS L 

dai 
d6 L 

dai 
9S L 



da i 



dS L 

dS N-l 

dS L 

d<t> 
dS L 
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This step was achieved in two stages: first, by adding each even row to the one above, and 
then by multiplying each even row by —1 and adding to it 1/2 the row above. 
Now note that 

d(p d(p da N dajq 



(3.9a) 

(3.9b) 

So we have 

•/ 



(3.10) 



dx 



da 



N- 



dciN dx 
da 



dx 



N-l 
2 



%a 7v 



- a-N 
1 

Oat 



3 I -77 



dx 



+ a 



dx 
daN 
dx 



la^ 



N-l 
2 



da i 


dai 




ddi 




da N _i 


dpi 




da,N-i 


daN-i 




dpi 


dBi 


da n+1 


da jy+i 



dpi 

daN 
dpi 



dBi 
daN 
ddi 



dai 
dp M 

dpM 



da n-1 
da n+1 

¥m 

daN 
dpM 



dai 

de M 
de M 



da N-l 
da n+1 

de M 

daN 
d0 M 



dai 
d5i 

d&i 



da N-l 

d&i 
da n+1 

dSi 
daN 
dSi 



dai 
d6 L 
daN-l 
d8 L 



da i 



dS L 
da n+1 



dS L 
daN 
d6 L 



So we have, where e is a quantity with modulus one that may vary at each occurrence, 



(3.11) 



e - 



N-l 
2 



dai 
dpi 



daN 
dpi 



dai 
ddi 



daN 
ddi 



dai 
dpM 



daN 
dpM 



dai 
96 M 



daN 
96 M 



dai 


dai 


dSi 


d6 L 


daN 


daN 


d&i 


d6 L 



For convenience let us set Uj 



(We thank the editor for suggesting how to improve 
the next part of the proof. Our argument was more lengthy.) The proof of the lemma is 
achieved by expressing the determinant ()3.11|) in terms of 



(3.12) 



dai 


dai 


dai 


dai 


dai 


dai 


d/3i 


dvi 


9Pm 


dv M 


dai 


da L 


daN 


daN 


daN 


daN 


daN 


daN 


dpi 


dvi 


9Pm 


dv M 


dai 


da L 



{-1) N A(/3i, z/i, . . . , p M , v m , oil)- 



To prove this equality note that the left hand side is a homogeneous polynomial of order 
N(N — l)/2 in the variables /3i, v\, . . . , [3m, ^m, Qti, . . . , oil- It is also antisymmetric under 
exchange of any two of the variables. Therefore, up to a constant, it must be proportional 
to the Vandermonde A(/3i, . . . , (3m, Vm, oci, . . . , «l) (defined in ()1.13j) ). The constant 
(—1)^ can be determined as follows. Let us denote by x\, . . . the variables z/ 1; . . . , 
(3m, v m, Oi\,...,OLL- The generic element in the Jacobian (|3.12|) is 

(3.13) H = (-l)V* 

where Cj^ is the j-th symmetric function in the xi excluding the variable x^. We want to 
compute the coefficient in front of 

(3.14) x^x^'-xn. 
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Now, any xj in any element of the Jacobian matrix appears at most with exponent one. 
Furthermore, for any given choice of j indices li, I2, ■ ■ ■ , lj, (l m 7^ k) the product 



(3.15) 



%ii ■£1% ' ' ' -El 



appears once and only once in the symmetric function e^fc. 
mial (|3.14j) can only come from the term 



It follows that the mono- 



(3.16) 



da 1 da,2 da^ 
dx\ dx 2 Oxj\[ 



in the expansion of the Jacobian and that the coefficient in front of it is (— i) Ar ( iV + 1 )/ 2 . The 
constant (— 1)^ is obtained by observing that the monomial (j3.14j) appear with a factor 
(_1)N(at-i)/2 j n eX p ans j on f the Vandermonde. 

The relation between the determinants ()3.11|) and ()3.12j) is obtained by analytically con- 
tinuing equations ()3.4|) in the complex planes of p and 9 and then by changing the variables 



from pi, 9x, . . . 
ther Jacobian 



Pm, m , Si,...,5 L to p 1 , ui, 



(3.17) 



dpi 

d0i 
.901 
d0i 



dpM 
d0i 

80m_ 

80i 
85i 

80i 



88 L 
80i 



dpi 
dvi 
&h 
dvi 



8pm 
dui 
89m 
dui 
881 
dui 



85 L 
dvi 



dpi 

80m 
ddi 

80 M 



dpM 
80u 
80 M 

80m 
8Si 

80 M 



8S L 
80a 



, Pm- 


V M , Oil, ■ ■ 


■ ,a L . 


dpi 


dpi 


dpi 


8vm 


dcei 


da L 


861 


ddi 


d6i 


8u M 


dai 


da L 


dpM 


dpM 


dpM 


8vm 


dcei 


da L 


89 M 




80m 


8vm 


dai 


da L 


881 


dSi 


dSi 


8vm 


dai 


da L 


8S L 


dS L 


dS L 


8vm 


dai 


da L 



Thus, we incur the fur- 



Using equations (|3.4j) the computation of this determinant is straightforward. The matrix 
is all zero except for 2 x 2 blocks on the diagonal for the first 2M columns, and then single 
elements on the diagonal for the last L columns. A given 2x2 block looks like 



(3.18) 



dpj 


dpj 




d0j 


dvi 




dOj 


de± 




d0 j 


dVj 






2u, 



and the single elements on the diagonal have modulus one. So, dividing ()3.12|) by ()3.17|) we 

have 



(3.19) 



Jr 



2 M~(N-l)/2 j TT J_ I 



A(/3i, fc, ■ ■ ■ , An, «i, . . . , a L ). 



Here e incorporates quantities that may depend on =-, /3 2 , =-,... , (3 m , =-, ai\, ...,«£ but 

01 02 Pm 

has absolute value one. This completes the proof in the case N is odd. 

When N is even, we want the Jacobian of the transformation from the real, independent 

variables Reai, Imai, . . . , Reajv , Imajv to the real independent variables p\, 9\, p 2 , O2, ■ ■ ■ , 

2 2 

Pm, &m, Si, ... , Sl- 
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So, we start with 



dReai 


9Reai 


dKeai 


dKeai 


9Rcai 


dKeai 


dpi 






90 M 


d5i 


9S L 


dlraai 


9Imai 


dlmai 


dlmai 


dlmai 


dlraai 


dpi 




dpM 


90 M 


d5i 


d5 L 


dKea jv 


dKea jv 


dKea jv 


dKea n 


dKea n 


dKea n 


~T 


T 


~T 




T 




dpi 


d8i 




m M 


d8i 


dS L 


dlraa jv 


dlraa jv 


dlraa jv 


dlraa jv 


dlraa n 


dlraa n 


T 


~T 




T 


T 


T 


dpi 


ddi 


dpM 


d0M 


d8i 


98 L 



(3.20) 



dai 
dpi 

dai 
dpi 



da jv 
X 



dpi 



da jv 
X 



dpi 



dai 



dai 
d0i 



da jv 

X 

96»i 



da jv 

X 

801 



dai 
dpM 

dai 
dpM 



da jv 
X 



c9a jv 
X 



dai 

90M 

dai 
d0 M 



da jv 
X 



9a jv 
X 



d<5i 

dai 
dSi 



da jv 
X 

05! 



da jv 
X 



da i 
95 L 

dai 
d6 L 



da jv 
X 



d<5 £ 



da jv 
X 



9p M 90 m d&i 

Now we want to transform the derivatives of Oj into derivatives of ajv-j, with the exception 
of ajv/2, which should give us derivatives of ajy. Note that has modulus 1 and 

CtiV/2 



Q-n/2 



a-N- 



da 



N 



(3.21) 

Therefore, we have 
(3.22) 

(3-23) _ 

a N/2 OX 

Thus we need to multiply the last row of ()3.20|) by 



dx 



d_ ( ON/2_\ 
dx \a N/2 ) 
1 da N/2 



O-N/2 



da 



N/2 



y^N/2 



dx 



-Ojv/2/(«Af/2) 2 and add to it l/ajv/2 times 
the row above. This procedure has multiplied the determinant by a factor — a^^/ (flN/2) 2 = 
— ajv/ajv/2- Note that a^ has modulus one, but ajv/2 does not. Thus, incorporating factors 
of modulus one into e, the Jacobian is 



(3.24) 



c — e ajv/2 



da 1 


dai 


9ai 


dai 


dai 


dai 


dpi 


80! 


dpM 


90 M 


dSi 


95 L 


c9ajv 


da jv 


c9<2 jv 


da jv 


dapf 


9ajv 


dpi 


90i 


dpM 


90M 


dSi 


9S L 



It is now just necessary to follow the steps as in the odd iV case in order to prove 
Lemma 13.11 □ 

4. Calculation of the Jacobian: real case 

We prove Lemma 12.21 It would be possible to cover a more general case, considering 
the number of zeros on the unit circle, the number in complex conjugate pairs located 
symmetrically with respect to the unit circle, and the number of real zeros. But our concern 
here is with the case that all zeros are on the unit circle so we only consider the case we require 
in this paper. We describe the odd degree case in detail and then discuss the modifications 
required for the even degree case. 
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We have a self-reciprocal polynomial f(z) that has a root at —1 and at /3%, /3 l , . . 
on the unit circle. Thus, f3j = The polynomial f(z) has order N = 2M + 1: 

f(z) = {z + l){z- Pi)(z -P 1 )---(z- f3 M )(z - P M ) 
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i P~Mi @m 



(4.1) 

We have the functional equation 
(4.2) 



z N + aiZ N-l + a2Z N-2 + 



f(z) = z N f 



a N _ 2 z + a N -iz + I. 



which implies the symmetry of the coefficients 

(4.3) a n = a N _ n . 

So the polynomial is determined by ai,. . . , au- We also have 

(4.4) a n = (-l) n e n (-l,Px,Pi, • • JmJm), 

where e n is the nth elementary symmetric function. 

We want the Jacobian of the transformation from the real variables a\,. . . ,a,M to the real 
variables t\ r . . ,tjy, where j3j 



e itj _ 



We start with 



(4.5) 



da± dai 


dai 


dti dt 2 


dt M 


da M da M 


daM 


dh dt 2 


dtM 



i M !3i---f3 



M 



dai 


dai 


dai 


dPi 


d(3 2 


d(3 M 


da M 


da M 


da M 


d/3i 


d(3 2 


d(5 M 



where we used the fact that 
(4.6) 



dctj da-; dfik 



dai 



dt k df3 k dt k d/3 k ' 
At this point our calculation becomes somewhat different than the complex case in the 
previous section. It is possible that similar methods could be used here, but the argument 
would be more sophisticated because (as we will see), the answer is not symmetric in all 
the variables. We will directly calculate the Jacobian, instead of deducing the form of the 
answer from its various symmetries. 

Following the method in for each m we have 



(4.7) 



-iy 



Pm 



"n— l.m 



4- p' 4- p' 



where we write e' nm for the nth symmetric function in —1 and all the /3j and (3^ = 

1, . . . , M, except for f3 m and with the convention that eo = e' m = 1 and e n = e 
if n is greater than the number of variables. Thus, we obtain 



so our determinant is 



da Ti 



1 



■'n—l.m 



M 



(4.9) 



n « 



771=1 



'0,1 
J 

"1,1 



"M-1,1 



"0,2 
5 1,2 



"M-1,2 



'0,M 
l l,M 



-M-XM 
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Here and following we use e to denote a number with absolute value 1, which may be different 
at each occurrence. 

Now we define the following polynomials: 

M „ N-2 



h{x) 



(4.10) 



[x 



M-l 



l)l[(x-(3 m )(x 



m= 1 



1 



E(- 1 ) n < 



X 



N-2-n 



n=0 



n=0 



M-l 



n=0 



where the next-to-last step used the fact that e' ni = — e' N _ 2 _ ne , since the roots come in 
conjugate pairs except for the extra root at —1. 
Note that fe(/3 m ) = if £ ^ m. Also we have 



M-l 



(4.11) 



A(e«) = 2S"-& £(-l)V B>€ cos((M - \ - n)t). 



n=0 



Therefore, with 
(4.12) V(t u ...,t M ) 





cos((M 




-cos((M 




( _ 1} m- 


- 2 cos(|*i) (- 


_1)M- 


_1 cos(|ii) 




cos((M 


- m 


-cos((M 


-i)*2) ••• 


( _ 1} m- 


- 2 cos(ft 2 ) (- 


_1)M- 


- 1 cos(|t 2 ) 




cos((M - 


- \)^m) 


-cos((M - 


- |)*m) ••• 


( _ 1} m- 


2 cos(|t M ) (" 


-1) M - 


1 OOs(|*Af ) 


we have 




























A (A) 












(4.13) 




V(t u . 


• , ^m)^r = 


Cm 





/2(A) 

























• • • Im(Pm) 







where 
(4.14) 



M / 

C M = e2" M ] [ i,„ - — 

m=l V ^ m 



Since 

(4.15) 2 sin(if) cos((M - | - n)t) = sin((M - n)t - sin((M - 1 - n)t), 
by multiplying each row by sin(|t„) we have that 

(4.16) V(t 1 ,...t M )-2 M Y[an(±t n ) 
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equals a determinant whose rows are 

(4.17) (sin(Mtj-) - sin((M - l)t,-) •■■ sin(2£j) - sin(i,) sin(t,)). 

By elementary column operations starting with the last column, this equals a determinant 
with rows 

(4.18) (sin(Mtj) sin((M - l)t,) ••■ sin(2^) sin(i,)). 
Since 

(4.19) sin((n + l)t) = sin(t) (2 ra cos n (t) + lower order terms in cos(t)) , 

by elementary column operations we find that the determinant whose rows are 1)4.18)1 equals 
a determinant with rows 

(4.20) sinfo) x (2 M - 1 cos M - 1 (t j ) 2 M ~ 2 cos^fo) ■•■ 2cos(t i ) l). 
Since this last determinant is a Vandermonde, we have shown 

V(t h . . . t M ) = 2 M ( M -^ TT JEW ACcosfa), . . • , cos(t M )) 

- LJ - sin i 

(4.21) " V2 ; 

= 2 M ( M " 1 )/ 2 J] cos(it n )A(cos(ti), . . . , cos(t M )). 

rt 

Combining lQ5|» . (l4"T4l . and (14"2T1) we have 

22) J R = £ 2~ M ( M + 1 )/ 2 . TimiP™ ~ Pm) Tim fm(Pm) 



n m cos(|t m )A(cos(*i), • • • , cos(t Af )) ' 



Since, 



m m k^m 

(4.23) = 2^~ M J] cos(|t m ) J] sin (^^) sin (hl±h 



k^£m 



and 



A(costi, . . . , cos^Af) = JJ(cosifc — COStj 



j<k 

(4.24) = 2 M ( M " 1 )/ 2 J] sin (^^j sin 

j<k 

and P m — P m = 2i sin(t m ), we find that, 

J K = e 2 M2 J] sin(t m ) J] sin {^) sin (**±^ 

(4.25) = e J](/3 m - J - ft) (ft - 

m j<fc 

as claimed. □ 
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In the case of even degree, only a few modifications are needed. The polynomials (|4.1()jl 
are replaced by 

m 1 



(4.26) 



m— 1 




x M-l-n + ^-M+l+n) + (-l)^ e ' M _ u , . 



SO 



M-2 

(4.27) /«(e«) = e ( 2 £(-l)V M cos((M - 1 - n)t) + (-1) M -V M _ L 

n=0 



The determinant V(ti, . . . , i/vf) has entries cos((M — 1 — n)tj), except for the last column, 
where n — M — 1, which is multiplied by a factor of |. Since 

(4.28) cos(nt) = 2 n_1 cos n (t) + lower order terms in cos(i), 

we recognize V(t±, . . . ,tju) as a Vandermonde, which is 2 M times smaller than the Vander- 
monde which appeared in the odd degree case. The only other differences in the calculation 
is to omit the factor 2 M Y\ sin(t n /2) from (j4.16|) and sin(tj) from (|4.2()|) . The overall effect 
of those factors is to multiply the Jacobian by 

(4.29) n#?%T = TTcos(^ 
v ; 1J -2sin% 11 \2 

That factor replaces riml^m + 1) = 2 M f] m cos(|t m ) which is now omitted from \[ m f m ((3 m ), 
the power of 2 making up for the factor of 2 M missing from V{t\, . . . , ^m)- So the end result 
is the exact same formula for the Jacobian. 

Appendix A. Epstein zeta functions associated to binary quadratic forms. 

An example ensemble of Dirichlet series with functional equation is the Epstein zeta- 
functions associated to a binary quadratic form. These functions have functional equations 
but no Euler product. Indeed it is conjectured (Sarnak's rigidity conjecture) that there does 
not exist a nontrivial continuous family of Dirichlet series with functional equation and Euler 
product. 

To create an Epstein zeta-function, begin with a binary quadratic form 
(Al) Q(m, n) = am 2 + bmn + en 2 . 

We assume that Q is positive-definite, meaning that Q(m,n) > if (m,n) ^ (0,0). This 
is equivalent to a > and A = b 2 — Aac < 0. In number theory one considers the case of 
a, b, c G Z, but we will allow a, b, c to be arbitrary real numbers. 
The Epstein zeta-function associated to Q is defined as 



< A2 > L «w : = E W 



Q(m, n) s 

where the notation means that we sum over all pairs of intgers (m, n) except for (0,0). 
This series converges for Re(s) > 1, and Lq(s) has a meromorphic continuation to the entire 
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complex plane except for a simple pole at s = 1. The Epstein zeta- function satisfies the 
functional equation 





(A3) ^ ns)L Q (s) = JLL-I r(l - s)L Q (l 



where T(s) is the Gamma function. Thus, Lq(s) has a functional equation similar to that 
of the Riemann (^-function and other L-functions of number theory. For certain quadratic 
forms the Epstein zeta-function can be expressed in terms of L-functions. For example, 

(A4) L m 2 +n2 (s) =4C(a)L(s,X4), 

where ((s) is the Riemann zeta-function and L(s, x&) is the Dirichlet L-function for the 
nontrivial character modulo 4. This Epstein zeta function is expected to satisfy the Riemann 
hypothesis, but note that its zeros are a superposition of two independent sets of zeros 
and so will not display the characteristic quadratic repulsion expected for L-functions with 
functional equation and Euler product. For the Epstein zeta function (jA4|) . and a handful of 
other special cases, Lq(s) has an Euler product. In the case that a,b,c are integers, Lq(s) 
is a finite linear combination of Dirichlet L-functions. However, for generic Q(m,n) there 
is no arithmetic structure to the Epstein zeta-function Lq(s), and hence no Euler product, 
and it is these "random Dirichlet series with functional equation" which illustrate that we 
are not studying an empty set of functions. 

There is a natural way to parametrize the Epstein zeta-functions. Write 



(A5) L Q {s) = a-s^T'- 



m + zn 



2s ■ 



where z = x + iy is a complex number satisfying \z\ 2 = c/a and 2x = b/a. Note that the 
conditions on a, b, c guarantee that such a number z exists, and z is uniquely determined if 
we also assume y > 0. The above formula for Lq allows us to relate the Epstein zeta-function 
to an Eisenstein series on SL(2, Z): 



(A6) L(z; S ):=V' 



y 



2^-^ \mz + n 



2s ' 



7)1,71 



Specifically, we have Lq(s) = 2(ay)~ s E(z; s). 

The Eisenstein series E(z;s) is a well-understood object from number theory. See ^2] 
for details. However, in that theory one usually fixes s and considers E(z; s) as a function 
of z. Indeed, as a function of z it has many fascinating properties: it is an eigenfunction of 
the hyperbolic Laplacian, and it is invariant (up to an automorphic factor) under the linear- 
fractional action of SL(2, Z) on the upper half-plane. Thus, the Epstein zeta-functions 
associated to a binary quadratic form are parametrized by the points in a fundamental 
domain for SL(2,Z): 

(A7) F=lz = x + iy : ~ < x < ~ y > 0, \z\ > 1 1 . 



In [5] evidence is given that such functions can be modeled by random trigonometric 
polynomials. 
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